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CN : Abstract 

^ ' A new treatment of the critical point of the two-dimensional interacting Bose gas is 

presented. In the lowest order approximation we obtain the critical temperature Tc ~ 
2Tm/[m\og{2Ti /mg)], where n is the density, m the mass, and g the coupling. This result 
is based on a new formulation of interacting gases at finite density and temperature which 

^ ■ 

■ is reminiscent of the thermodynamic Bethe ansatz in one dimension. In this formalism, 

the basic thermodynamic quantities are expressed in terms of a pseudo-energy. Consistent 
resummation of 2-body scattering leads to an integral equation for the pseudo-energy with 

"> 

• a kernel based on the logarithm of the exact 2-body S-matrix. 
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I. INTRODUCTION 



The properties of interacting Bose gases can be very different depending on the 
spatial dimension. This has become especially interesting in recent years due to 
the possibility of experimentally realizing lower dimensional cold gases with mag- 
netic and optical traps. In 3 dimensions there is a critical point for Bose-Einstein 
condensation (BEC) even in the non-interacting theory, with a critical temperature 
= 47r (n/C(3/2))^/^ where n is the density and ( Riemann's zeta function. In two 
dimensions the same formula becomes Tc = 47™/C(l) and since ({1) diverges there 
is no critical point at finite temperature for the non-interacting gas. It is believed 
however that the two-dimensional interacting gas has a critical point in the univer- 
sality class of the Kosterlitz-Thouless transition rather than BEc|l|, 2], [s], 
review see jsl.) This transition has recently been observed in experiments 

In the one-dimensional case the particles effectively behave as fermions, the model 
is integrable, and an exact solution is known based on the thermodynamic Bethe 



41. (For a 
3. 



ansatz (TBA) 



Sj. The one-dimensional case illustrates the importance of non- 
perturbative methods like the TBA for understanding the physics. Whereas the 
usual perturbative finite-temperature Feynman diagram method involving Matsub- 
ara frequencies in loops is a standard approach which entangles zero temperature 
perturbation theory with quantum statistical sums, the TBA represents an entirely 
different organization of the free energy which essentially disentangles the two. More 
specifically, the only property of the model it is based on is the exact 2-body scat- 
tering matrix computed to all orders in perturbation theory at zero temperature. In 
principle a TBA-like organization of the free energy is possible for non-integrable 
theories in any dimension based on the formula in 9| which expresses the partition 
function in terms of the exact S-matrix. The derivation of the TBA as given by 
Yang and Yang however was not based on the result in 9|] but rather relied on the 
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factorizability of the multi-particle S-matrix into 2-body S-matrices for integrable 

systems. For non-integrable systems the formalism in 9| can be quite complicated 

since the free energy contains N-body terms which do not factorize. Nevertheless, for 

a dilute gas the consistent resummation of the 2-body scattering terms can represent 

a useful approximation that is intrinsically different from other methods. 

In this work we derive explicit expressions for the free-energy, occupation numbers, 

etc., in this 2-body approximation for non-integrable models in any spatial dimension. 

The result is TBA-like: everything is expressed in terms of a pseudo-energy that 

satisfies an integral equation with a momentum-dependent kernel which is related 

to a matrix element of the logarithm of the S-matrix. Our analysis builds on the 
□ 

previous work[10l|. and the final result presented here contains several important 
technical improvements. 

Most of this paper is devoted to developing the formalism in generality. In section 
III we describe contributions to the free energy with diagrams, not to be confused 
with finite temperature Feynman diagrams, where the vertices represent A^-body 
interactions for any A^. The cluster decomposition of the S-matrix is necessary to 
establish the extensivity of the free energy. In section IV we derive an integral 
equation for the occupation numbers (filling fractions) from a variational principle 
based on a Legendre transformation that exchanges the chemical potential with the 
filling fraction. In section V we present the integral equation that consistently resums 
the infinite number of 2-body diagrams. In section VI the 2-body kernel that appears 
in the integral equation is derived for interacting Bose gases in 2 and 3 dimensions. 
In section VII we compare our approximation to the exact TEA for the 1-dimensional 
case. 

This new formalism is illustrated in the 2-dimensional case. Since the 2-body 
approximation is reasonably simple in its final form, we present a self-contained 
analysis of the critical point of the Bose gas in the next section, where we derive an 



expression for the coupling constant dependent critical density. 



II. CRITICAL POINT OF THE TWO-DIMENSIONAL BOSE GAS 

In this section we illustrate the main ingredients of our formalism by applying it 
to the critical point of the two-dimensional Bose gas. Our treatment is significantly 
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different from previous onespj, 12, 13|, U]. Although we find results that are similar to 
the known results, we find our treatment to be considerably simpler and transparent 
and doesn't rely on an effective description of vortices. As will become clear, the 
full power of our formalism is actually unnecessary to obtain these results when the 
coupling is very small. We set h = ks = 1 unless otherwise indicated. 

The interacting Bose gas in 2 spatial dimensions is defined by the following hamil- 
tonian: 

H = J d^x(^^|W(x)r + ||<^(a;)r) (1) 

where is a complex field and m is the mass of the particles. In two dimensions, the 
combination mg/h'^ is dimensionless; in the experiment {g] it is approximately 0.13. 

In the formalism developed in the sequel, the occupation numbers, or filling frac- 
tions, are expressed in terms of a pseudo-energy e(k), and the density has the fol- 
lowing form: 

J (27r)2 e/^^W - 1 ^ ' 

where (5 is the inverse temperature T. In an approximation that consistently resums 
all the 2-body interactions, the pseudo-energy satisfies the integral equation in eq. 
(l66l) where the kernel 6*2 (k, k') is related to a logarithm of the exact 2-body scattering 
matrix and given explicitly in eq. flHOl) . for convenience reproduced here: 

„ , , 8 / mg/8 , 
G2(ki,k2) = arctan 7 ^ (3) 

Vl + i^log(^ 



where A is an ultra-violet cut-off. When the coupling g is small, the kernel is small 
and the eq. fl66|l can be approximated by: 

where /i is the chemical potential and cuk = k^/2m is the single particle energy. 

Since when (7 = 0, e(k) = c^k — /i, an approximate solution to (j4]) is obtained by 
substituting e(k) = uj\f^ — ^ on the right hand side. As in Bose- Einstein condensation, 
the critical point is defined to occur at the value of the chemical potential /ic where 
the occupation number at zero momentum diverges, i.e. e(k = 0) = 0. This gives 
the following integral equation for 

At very weak coupling, to leading order we can neglect the momentum dependence 
of the kernel and simply take G2 = —Q- Consider first the case where the interaction 
is attractive, i.e. g The equation ([5]) becomes 

/5/ic = -^log(l-e^''^) (6) 

For negative g, the above equation has negative /ic solutions, which makes physical 
sense since energy is released when a particle is added. Here the attractive interaction 
makes a critical point possible at finite temperature and density. The critical density 

Suppose now the attraction is repulsive, i.e. g > The equation ^ now has no 
real solutions for ^c- This can be traced to the infra-red divergence of the integral 
when changes sign. Let us therefore introduce a low-momentum cut-off /cq, i.e. 
restrict the integral to |k| > k^. The equation ([5]) now becomes 

= ^ [/3eo - P^^c - log (e^^^""^^) - l)] (7) 



where eo = kl/2m. The above equation now has solutions at positive /Zc, however for 
arbitary there are in general two solutions. If is chosen appropriately to equal a 
critical value kc-, then there is a unique solution /ic. This is shown in Figure [Tj where 
the right and left hand sides of equation (JTj) are plotted against /ic for two values of 
/co, one being the critical value kc- Since the slopes of the right and left hand sides 
of eq. ([7]) are equal at the critical point, taking the derivative with respect to /ic of 
both sides gives the additional equation: 

/3(e, - /ij = log (l + ^) (8) 

where ec = k1/2m. The two equations ([71 [8]) determine /ic and ec- 

+ (9) 
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Note that €c goes to zero as g goes to zero. Finally, making the approximation 
e(k) ^ Ld^ — fic in eq. ([2]) and using eq. ([7]) one obtains the critical density Uc ~ yUc/fi': 



nc~^-|^log 1 + — (11) 
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The above can be expressed as ricA^ = log(l + 271 /mg) where A = h\j2TT /mksT is 
the thermal wavelength. As expected, at fixed density, the critical density goes to 
infinity as g goes to zero. 

In obtaining the above result we have only kept the leading term in the kernel G2 
at small coupling, which is momentum- independent. Systematic corrections to the 
above results may be obtained by keeping higher order terms in the kernel. However 
since the kernel involves logk^, at small coupling these corrections are quite small, 
as we have verified numerically. 

In numerical simulations of the critical point, it was found that /ic ~ 
log(2C/m(/) where C is a constant^. This form agrees with our result ([9]) 




FIG. 1: Left and right hand sides of eq. d?]) as a function of (3fic at two values of the 
infra-red cut-off eo, one being the critical value €c- 

in the limit of small g, however, whereas we find C = ir, numerically one finds a 
significantly higher value C ~ 13.2. Assuming that C was not overestimated numer- 
ically, one possible explanation for this discrepancy is based on the renormalization 
group. The beta-function for g can be found by requiring that the kernel G2 in eq. 
(l80l) be independent of the ultra-violet cut-off A. To lowest order this gives: 

= ^ (12) 
dA 4:71 ^ ^ 

which implies the coupling decreases at low energies. Integrating this equation leads 

to the cut-off dependent coupling g{A): 

2tx 271 ( mgo 



(l + ^log(Ao/A)) (13) 



mg{A) mgo 

where go = 5'(Ao). Thus at low energies where A < Aq, the combination 27i/mg 
grows. However since the dependence on the cut-off is only logarithmic, it is not 
clear that this could explain the large C discrepancy. This seems to imply that 
higher N-body interactions are important. 



III. PARTITION FUNCTION IN TERMS OF THE S-MATRIX 



A. Formal expression for Z in terms of the S-matrix 

We are interested in the partition function: 

= Tre-^(^-'^^) (14) 

where /3 = l/ksT and fi is the chemical potentiaL We henceforth set ks = 1. The 
trace will be performed over the free particle Fock space. Assuming only one kind of 
particle, the trace is computed based on the following resolution of the identity: 

oo 

1 = 5^^ /(ciki)(rfk2)---(dkjv) |ki..-kjv)(ki..-k^| (15) 

where (dk) = (i''k/(27r)'^ and d is the spatial dimension. We adopt the following 
conventions: 

|ki,...,k^) = ai^---4^|0) (16) 



ak4, -s ai,ak = (27r)'^5("')(k-k') (17) 

where s = 1 corresponds to bosons and s = —1 fermions. For simplicity we will use 
the following notation: (5k,k' = ^'^'^''(k — k'). The energy of a free one-particle state 
will be denoted as cuk- We will also need: 

= 5(^n0) = (18) 

where V is the volume. 

As usual, we assume the hamiltonian can be split into a free part Hq and an 
interacting part: H = Hq + Hj. Let S{E) denote the S-matrix operator in the formal 
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theory of scattering where E is an off-shell energy variable. It can be expressed in 
the conventional manner: 

S{E) = l + 2m6{E-Ho)f{E) (19) 

where Hq is the free hamiltonian operator and everywhere the over-hat de- 
notes a quantum operator. On shell the matrix elements of T[E) are T = 
(27r)'^(5{k,k'}-^{k}^{k'} where A4 are the scattering amplitudes. 
Our starting point is the basic formula derived in 91]: 

Z = Zo + ^ (iEe-^^Tr (s^^ HeS) (20) 
where XdsY = X{dEY) — (dF;X)Y, and Zo is the free partition function. We will 



work with the equivalent expression 



Z = Zo + —/ rf^e"^-^ Trim 9e log S(E) (21) 

271 J 

Since log S'(i?) oc 6{E — Hq), let us define the operator W as follows: 

Im log S{E) = 27r6{E - Ho) W{E) (22) 

Integrating by parts, one then has 

Z = y"rfEe-^^Tr(5(E-Ho)(W(E)/? + l)) (23) 

The advantage of the above expression is that in taking the trace, Hq can be replaced 
with the free particle energies and the integral over E performed. For simplicity of 
notation, define 

W = W(3 + 1 (24) 
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B. Cluster decomposition 



As we will see, clustering properties of the S-matrix ensure that the free energy 
depends only on connected matrix elements and is proportional to the volume. For 
any operator X, the cluster decomposition can be expressed as follows: 

(*'|X|^)= {%\X\^i)c{%\X\^2)c{%\X\^3)c--- (25) 

partitions 

where the sum is over partitions of the state into clusters |^'2)7 (The 

number of particles in and is not necessarily the same.) The above formula 
essentially defines what is meant by the connected matrix elements {X)c- 

In order for the free energy to be extensive, W must cluster properly, as we now 
describe. Introducing the short-hand notation |12 ■ • ■ ) = |kik2 ■ • ■ ), for the operator 
W the above equation reads 

= (26) 
{1'2'\W\12) = {1'2'\W\12), + {l'\W\l)c{2'\W\2), + s{2'\W\l)c{l'\W\2)c 
(l'2'3'|i?|123) = {l'2'3'\W\123)c 

[l'\W\l)c{2'3'\W\23)c + permutations 

{l'\W\l)e{2'\W\2)c{3'\W\3)c + permutations 

etc. The subscript [*]9 denotes the number of terms within the bracket, so that for 3 
particles there is a total of 1 + 9 + 6 = 16 terms on the right hand side of the above 
equation. The connected matrix elements are characterized by being proportional to 
a single overall (5-function. 

When {k'} = {k}, the above cluster expansion leads to another specialized type 
of cluster expansion that is suited to the computation of the partition function. From 
eq. ( !26l) we define ■u;Ar(ki, ...,k7v) as factors that cannot be written as a product of 
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separate functions of disjoint subsets of the kj: 



{1\W\1) 
(12|W^|12) 
(123|W^|123) 



(l|W^|l), = iI;i(ki) 
W2{ki, ka) + ^;i(ki)il;i(k2) 



fi;3(ki,k2,k3 



W2(k2, k3)u}i(ki) + W2(ki, k3)il;i(k2) 
+ M;i(ki)wi(k2)wi(k3) 



W2('kl,'k2)Wi 



etc. The above definition imphes that the Wn are sums of terms involving the 
connected matrix elements of W. For instance: 



w;2(ki,k2) = (12|W^|12), + s(2|W^|l),(l|iy|2), 
il;3(ki,k2,k3) = {123\W\123)c+\s{2\W\l)c{13\W\23)c + verm. 



(28) 



+ 



{3\W\l),{l\W\2),{2\W\3), + peTm. 
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A combinatoric argument then shows: 



N=l ' 



(29) 



where z = e^^ and oji = tUk^. 

All of the w are proportional to an overall {2ttY6^'^'^ (0) = V, so let us factor out 
the volume V and define: 



WN{ki, ...,kjv) = VwN{ki, ...,k 



N) 

The free energy density is then completely determined by the functions wn'- 
^ = -— log Z = -- — y (rfki) ■ ■ ■ (dk^) e-^^«-» w^{k,, k;v) 



(30) 



(31) 



For free particles, the only non-zero connected matrix element is the 1-particle to 
1-particle one: (2|iy|l) = (2|1). This still implies Wn is non-zero for any N: 



(uNil, 2, .., N) = s^'-'iN - 1)! (N|l)(l|2)(2|3) ■ ■ ■ (N - 1|N) 



(32) 
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Note that although wn is a product of N 5-functions, it only has a single 5(0). All but 
one k integral is saturated by 5-functions at each N and one obtains the well-known 
free-particle result: 

C. Diagrammatic description 

In the interacting case, terms involving the one-particle factors (2|W^|l)c can be 
easily summed over as in the free particle case. We assume that |k) is a stable 
one-particle eigenstate of H with energy uj^. 

(k'|l?|k), = (k'|k) = (27r)'^<5k,k' (34) 

Consider for example the W3 contribution coming from the second terms in eq. fl28l) : 

logZ = .... + ^ / ((iki)(rfk2)(c/k3)e-^^-i"» (s(2|l?|l)e(13|l?|23), + perm.) + ... 

(35) 

Because of the 5-function, one can do one k integral. Combined with the primary 
term coming from W2 one finds: 

log Z = ... + — / (dki)(dk2)e-^(^^+'^2) (1 + S2(e-^'^i + e-^'"^)) (12|W^|12), + ... 

(36) 

It is easy to see that the terms in the higher wn that involve one ■W2 and — 1 wis 
contribute additional terms in the parentheses of the above equation that sum up to 
]^^^-,^(l — sze~'^'^')~^. To summarize this result in the most convenient way, define 
Y^'j^ as the sum over Wn terms minus any terms involving the 1-particle factors Wi. 
Then ^ 

•^ = -^o-^E'^ / {fl{dK)fo{kA w^ik„...,k^) (37) 

^ N>2 \i=l / 
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where JFq is the free contribution in eq. (l33|l and we have defined the free "filhng 
fraction" : 

/"(k) - (38) 

The functions W]s[ have many contributions built out of the connected matrix 
elements of W . Define the vertex function V(k'^, k^; ki, k„) as follows: 

(k;, Kn\W\ki^ K)c = /5(2vr)'^4,k' V(k;, \^^- ki, K) (39) 

The vertices V have no temperature dependence and are essentially on-shell matrix 
elements of log S. To lowest order in the scattering matrix T, this vertex is just the 
scattering amplitude: V = M.. In the sequel we will compute V2 for the interacting 
Bose gas. Given our present interest in non-relativistic systems, only the matrix 
elements with n = m are of concern here. We highlight this fact by attaching a 
subscript n to V: 

Vn{K k:,; ki, k„) = V(k;, k^; ki, K) (40) 

Each term in the corrections T — Tq can be represented graphically. We stress 
that this diagramatic expansion as defined here, while being analogous to the usual 
Feynman diagrams, is completely unrelated to the (finite temperature) Feynman 
diagram approach. This is clear since each vertex in principle represents an exact 
zero-temperature quantity summed to all orders in perturbation theory. The correc- 
tion JF — jFg is the sum of all connected "vacuum diagrams" , made of oriented lines 
and vertices with n incoming and n outgoing lines attached, n being any positive 
integer. Note that there must be at least one vertex; i.e. the circle isn't allowed since 
it is already included as jFg. A diagram is then evaluated according to the following 
rules and represents a contribution to T: 

1. Each line carries some momentum k. To such a line is associated a factor of 
/o(k). 
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2. To each vertex, with the incoming set and the outgoing set of momenta be- 
ing {ki, . . . , k„} and {k'^, . . . , k^} respectively, associate a factor of the vertex 
function {27ry(3 V„(k;, k^; ki, k„). 



3. Enforce momentum conservation at each vertex. 

4. Integrate over all unconstrainted momenta with o?'^k/(27r)'^. 

5. Divide by the symmetry factor of the diagram, defined as the number of per- 
mutations of the internal lines that do not change the topology of the graph, 
including relative positions. This factor is identical to the usual symmetry 
factor of a Feynman diagram. 

6. For fermions, include one statistical factor s = — 1 for each loop. 

7. Divide the result by the volume of space V = {2nY6{0) and by /5. Note 
that this is equivalent to dividing the space-time volume V/S since at finite 
temperature, (3 is the circumference of compactified time. 

The structure of this diagrammatic expansion is very similar to that in the work of 



Lee and Yang[ll|. however the vertices are different in the two approaches. 

To illustrate this diagrammatic expansion, we will now describe the first few terms. 
There is only one diagram that has two lines, as shown in Figure [2l This diagram 
corresponds to: 

^= ....^ J (dki)(rfk2)/o(ki)/o(k2)V2(ki,k2;ki,k2). (41) 

Comparing with the expression fl37|) . we identify W2(ki, k2) = /3V2(ki, k2; ki, k2). At 
the next order, there is again only one diagram consisting of three lines, shown in 
Figure El and this implies W3(ki, k2, ks) = s/5V3(ki, k2, ka; ki, k2, ks). 
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(a) 




FIG. 2: (a) Contribution to W2- (b) Contribution to ^3. 






FIG. 3: Diagrams contributing to W4, 

There are multiple diagrams with four lines, one of them includes V4, and the 
others can be built from V2. They are shown in Figure [3l Their sum is 

4 



T=...+ 



n(^k,)/o(k, 



i=l 



— V4(ki,...k4;ki,...k4) 



(42) 



+ -{2TlY5i\^2 - k3)/3V2(ki, k2; ki, k3)V2(k2, k4; kg, k4) 

4 

+ ^(27r)'^5( V k,)/52V2(ki, k2; k3, k4)V2(k3, k4; ki, k2)" 

From this expression, we can identify 1^4, which is now a sum of terms built from 
the vertices V4, V2. 

An interesting class of diagrams are the "ring-diagrams" shown in Figure 3, since 
they can be summed up in closed form. Let us define 



G'2(k,k')=V2(k,k';k,k') 



(43) 
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n 



FIG. 4: Ring diagrams. 

A ring diagram with external rings attached will have the following value according 
to the diagrammatic rules: 

1 r ^ 

— j (rfko) \{{d\^i) [/5G2(ko, k,)/o(k,)] /o(ko)^ if TV > 1; (44) 



^ j (rfko)(dki)/3G2(ko,ki)/o(ki)/o(ko) if Ar = 1. 



Apart from the A^ = 1 term which has a different coefficient, all other terms can be 
identified with the power series for the logarithmic function. This class of diagrams 
are then resummed into: 



•J ring 



- y (rfko) log [\ - s/o(ko) J (rfki)/3G2(ko, ki)/o(ki) j (45) 



^ / (dko)(rfki)G2(ko,ki)/o(ki)/o(ko) 
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IV. LEGENDRE TRANSFORMATION AND INTEGRAL EQUATION 



A. Generalities 

Given one can compute the thermally averaged number density n: 

n = -^^l (dk) /(k) (46) 

The dimensionless quantities / are called the filling fractions or occupation numbers. 
One can express JF as a functional of / with a Legendre transformation. Define 

G = J^{lj)+fin (47) 

Treating / and /i as independent variables, then using eq. (146!) one has that d^G = 
which implies it can be expressed only in terms of / and satisfies 6G/6f = fi. 
Inverting the above construction shows that there exists a functional FifyfJ^) 



r{f,fi) = G{f)-f, J{dk)f{k) (48) 



which satisfies eq. (l46ll and is a stationary point with respect to /: 

1=0 (49) 

The above stationary condition is to be viewed as determining / as a function of fi. 
The physical free energy is then = F evaluated at the solution / to the above 
equation. We will refer to eq. (149!) as the saddle point equation since it is suggestive 
of a saddle point approximation to a functional integral. The function F is also 
required to satisfy 

|^ = -/w/(k) (50) 

In the sequel, it will be convenient to trade the chemical potential fi for the 
variable /q. We will need: 

^ = PMl + sfo) (51) 
17 



In the free theory the functional Fo is fixed only up to the saddle point equations, 
thus its explicit expression is not unique. As we will see, the appropriate choice that 
is consistent with the diagramatic expansion is: 

ro(/, /o) = / (^k) (^s log(l + ./) - (/^) ) (52) 

(This is different than the choice made in [l^.) The saddle point equation = 
implies / = /o, and inserting this back into F o gives the correct free energy JFq for the 
free theory. Furthermore, F o satisfies eq. f l50|l if one uses the saddle point equation 

f = U 

Let us now include interactions by defining 

F = Fo + Fi (53) 

where To is given in eq. (!52|) and we define U as the "potential" which depends on 
/ and incorporates interactions: 

ri = ~l (dk) Uifik)) (54) 

B. Integral equation 

The physical filling fraction / can be obtained by differentiating JF in (!37|) with 
respect to the chemical potential /i, and the result is 

" N 



oo „ 

/(k) = /o(k) + /o(k) (./o(k) + 

N=l ' 



i=l 



ty7v+i(k,ki, . . . ,k7v). 

(55) 

Since /o is represented by a line in our graphical expansion, it may be tempting to 
postulate that / would be a fully dressed line, but this interpretation isn't exactly 
correct. Graphically, a dressed line is obtained by summing over all possible inser- 
tions due to interactions. The sum can be generated in the following way: take a 
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(connected) vacuum diagram, cut an internal line into two external legs, and the 
resulting diagram contributes to the sum. If we sum over all possible ways of cutting 
an internal line from a vacuum diagram, and furthermore sum over all such diagrams, 
we have generated all correction terms to the so-called propagator. Cutting a line 
and replacing it with two legs corresponds to replacing a factor of /o(k) by s/o(k)^. 
We have the extra factor of s because cutting one line removes one loop integral. 

Recall that JF — JFq is the sum of all connected vacuum diagrams. The action of 
the derivative d/dfi makes sure that we've picked out every internal line in every 
diagram. However, it replaces /o by /o(l + s/o). Consequently, we define the dressed 
line 7(k) by 

7(k) = /o(k) + (-4^^) (/(k) - /o(k)) . (56) 



,s/o(k) + l. 

It will be useful to define ek as the single particle pseudo-energy at momentum k as 
follows: 



1 

Then the eq. (!56|) can be equivalently expressed in the simpler form: 



/(k) = e/'(^k-<^k+/^)y(k) (58) 
The following identity will be useful: 

||il = (2.)^*(k-k0-^^ (59) 

5/(k') 1 + s/o(k) 

This fully dressed line should satisfy an equation analogous to the Dyson equation 
for the full propagator: 

7(k) = /o(k) + ./o(k)/(k)E(k), (60) 

where S(k) is the sum of all (amputated) one-particle (IPI) insertions. Figure [5] 
represents this equation graphically. 
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FIG. 5: Graphical representation of the integral equation (|60p . The double line represents 
/, and the single line represents /q. 

Again we consider the analogy to Feynman diagrams. Take a two-particle ir- 
reducible (2P1) vacuum diagram with each internal line being fully dressed /, and 
remove one internal line to form an amputated IPl insertion. Summing over all ways 
of removing a line of a diagram, and then summing over all 2P1 vacuum diagrams, 
the result is S. 

Let — /5F 1 be the sum of all 2P1 vacuum diagrams, with internal lines being / 
instead of /q. The operations above can be neatly summarized by the following: 



The factor of s is inserted by hand because the removal of one internal line decreases 
the loop count by one. One correctly surmises that F i turns out to be the interaction 
part of F- Indeed, when there is no interaction, all vertex functions are zero and F i 
vanishes. The "Dyson" equation now reads: 



and is actually an integral equation for /. 

If we demand that the saddle point equation of the functional F leads to flB21) . 
then F o can be uniquely fixed up to an additive constant and an overall scale. If we 
further demand that Fo reproduces the result of the free theory, then one can show 
that the only consistent choice is F = Fo + Fi, where To is given in eq. fl52|) and 




(61) 



/(k)-/o(k)+/3/(k)(l + ./o(k)) 



SFi 



0. 



(62) 



5/(k) 



ri = ^^^^(/o - /) 



(63) 
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i.e. is given by the 2PI diagrams with /q replaced by the fully dressed /. In 
Appendix A, we prove that the saddle point value of F is exactly JF, thus completing 
the derivation. 

V. TWO-BODY APROXIMATION 

We shall now make the following approximations: 

1. include only the four-vertex, i.e. V2, in our diagrammatic expansion; 

2. consider only the foam diagrams, shown in Figure 

This should be the leading contribution in the low-density, high-temperature limit. 
When the density of the particles is low, we expect the two-particle scattering to be 
the most significant. At any given order of /o, a foam diagram has the least number 
of vertices; since each vertex contributes a factor of (3, the foam diagrams dominate 
in the limit of small (3. In practice, this approximation amounts to setting 



and truncating all terms with more /. 

Within this approximation, the integral equation fl^2]) can be rewritten into a 
self-consistent equation for the single-particle pseudo-energy. With the expressions 
(l64l) and (!59l) . the equation (162|1 reduces to 



In terms of the pseudo-energy ek defined in eq. (j57l) , using eq. (l58l) one can rewrite 
(l65l) as follows: 




(64) 




(65) 




(66) 
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FIG. 6: Foam diagrams. 



One may also want to find the free energy density of tlie system in terms of the 
physical filling fraction /. By using (165!) and (!56|) . it can be shown that 

ri = -\j (rfk)(dk07(k)7(kOG2(k,k') (67) 

-2(3] ^^^^ 1 + ./o(k) • 
The free energy density then takes the simple form: 

T = Fo + Fi (68) 
^ ^(rfk)fslog(l + s/) ^ ( f-f' 



In summary, the two-body approximation consists of the integral equation (166|) 
for the pseudo-energy defined in eq. ( 1571) . and the above expression ( l68l) for the 
free energy. The kernel that enters this approximation will be computed in the next 
section for Bose gases. 
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VI. TWO-BODY KERNEL FOR INTERACTING BOSE GASES 



A. The kernel in arbitary dimensions. 

The interacting Bose gas in d spatial dimensions is defined by the hamiltonian 
H = jd'^ (^|V0(x)r + (69) 

In this section, we derive the two-body kernel G'2(ki)k2) that enters the approx- 
imation of the last section. According to the definition in Section III, the kernel is 
the connected part of the diagonal matrix element of ImlogS'(£'). Since unitarity of 
S implies it is a pure phase, one has 

2ti5 {e - ^(k? + k^)^ V G2(ki,k2) = -^(ki,k2| log^(E)|ki,k2). (70) 

where V is the spatial volume. To define log(5'), we first follow the usual convention 
and write the two-body S-matrix as the sum of the identity operator and the T- 
matrix, as in eq. (fT9|) . The matrix elements of T are the following: 

(pi,P2|f (E)|qi,q2) = (27r)'^5(pi + p2 - qi - q2)A?,„,,^p„p,(E). (71) 

where M.q^^q^^p^^p^{E) denotes the two-body scattering amplitude from the in-state 
|qi, q2) to the out-state |pi, P2), assuming the energy E of the in-state takes on some 
off-shell value. We'll denote the on-shell version by M.q^^q^^p^,p^. Now log(S') may 
be defined as a series of T: 

\og(S{E)) = -y2 ^ f27r(5(E - Ho) i f{E)\ " . (72) 

n=l 

The connectedness is guaranteed because, in terms of Feynman diagrams, T does 
not contain any disconnected pieces. 
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The diagonal matrix elements of the n-th term of the above operator series, 



(ki,k2 



2ti5{E-Hq)iT{E) |ki,k2) 



can be found by inserting copies of the resolution of the identity. 

Note that the leftmost energy 5-function hits the state (ki, k2|, and consequently 
puts the energy E on-shell. Moreover, given that the S-matrix is invariant under 
rotations and Galilean boosts, the on-shell amplitude M. can only depend on the 
difference of the incoming momenta, which remains invariant due to conservation of 
energy and momentum. The amplitude M. therefore remains constant throughout 
the entire available phase space, and can be taken out of any integration over the 
phase space. All put together, the diagonal matrix element of the n-th term of the 
series ( 1721 ) is 

1 



(ki,k2 



2n6{E - Ho) I T{E) |ki, k2) ={27r)6 { E - —{kf + k',) ] V 



(73) 



x(zM(|ki-k2|))"J"-\ 

Here X = X(|ki— k2|) is the volume of phase space available to the scattering out-state 
resulting from summing over all two-particle states in every intermediate resolution 
of identity: 

J(|ki-k2|) = ^ |(rfp^)(rfp,)(2vr)'^+i5(^i^__L(p2 + p2)^5W(p_p^_p^)^ (74) 

where E = (k^ -|-k2)/2m and P = ki + k2 are the total energy and momentum of the 
in-state. Using this in (172|) . the resummation goes through trivially, and one obtains: 

(ki, k2| log(^(E))|ki, k2) = 27r5 (^E - ^(k? + k^)^ ^ log (1 + tIM{\k, -k,\)) 

(75) 

Finally, comparing this with (ITOl) . we have arrived at 

G'2(ki, ks) = ^ log (1 + zIM{\k,-k2\)). (76) 
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FIG. 7: Ladder diagram resummation of the two-body scattering amplitude. 

The exact ampUtude Ai in any number of dimensions d can be calculated by 
resumming ladder diagrams (Fig [7]): 



zM = {-ig) + \{-igf j {\oov) + \{-igf (^j (loop)^ 



-tg[l + '-l I (loop) ' 



(77) 



where the loop integral is 



(loop)= Iid-)id^){,_j,y2^^^^) [j^^Z^y^i^^^y^^e^ 



B. The kernel in two dimensions. 

The first step of the calculation is to compute the kernel by setting d = 2 in 
equations (j76l) . ( j74l) and ( jTTl) . In dimension d = 2 the phase space volume X = m/4. 

The loop integral diverges ind = 2, and must be regularized with a UV momentum 
cutoff A. The final form of is: 

M{\k,-k,\) = \ ^ (79) 

When all put together, the kernel G2 is the following: 

lki-k2 



G2(ki,k2) = log 



mg/S 



(80) 



— arctan 
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C. The kernel in three dimensions 



Computation of the two-body kernel in ci = 3 has httle different from the 2- 
dimensional case. The phase space volume fl71|l turns out to be: 

X=— |ki-k2|. (81) 

OTT 

Once again the loop integral is UV divergent and in principle needs a cutoff A. The 
divergence is linear in A, and can be completely absorbed by a redefinition of the 
coupling constant. We proceed to define the renormalized coupling by 

- = i + (82) 

9r ^ 47r2 

The amplitude (1771) becomes 



The two-body kernel fl76l) is then 



G'2(ki,k2) = -i I, , I log 
m ki — k2 



1 - 


imgn 
Wn 


|ki 


-k2| 


1 + 


imga 
167r 


|ki 


-k2| 



16vr f!!^^|ki 



(84) 



arctan 



mlki — k2| V 167r 



VII. ONE-DIMENSIONAL CASE: COMPARISON WITH THE THERMO- 
DYNAMIC BETHE ANSATZ. 

In one spatial dimension, the thermodynamic Bethe ansatz (TBA) is applicable 
if the model is integrable, leading to an exact expression for the free energy. As 
explained in the Introduction, our formalism is modeled after the TBA, and is con- 
structed from the same ingredients: the filling fractions and free energy are expressed 
in terms of a pseudo-energy e(/c), and the latter satisfies an integral equation based 
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on the exact S-matrix. If the model is integrable, the n-body S-matrix factorizes 
into 2-body S-matrices, and the TBA incorporates these n-body interactions in its 
derivation and the final result is expressed only in terms of the 2-body S-matrix. The 
TBA thus goes beyond the 2-body approximation described in section V. Neverthe- 
less, our formalism should coincide with the TBA to lowest order in the scattering 
kernel, and in this section we show this is indeed the case. 



A. Kernel 



In this section, we shall specialize in the case d = 1. Also, we shall set m = 1/2 
throughout this section. The on-shell value of energy for the two-particle in-state 
\ki, is E = kl + kl- With this identification, the phase space volume fl7^ turns 
out to be 

'''' 

The two-body amplitude (1771) in d = 1 is 

= ~\ . (86) 

+ 4(fci-fc2) 

The kernel (176ll then becomes: 

G.(.....)^-2.(..-..).og( ;^;:g;;^/^ ). (87) 



B. Comparision with TBA 

In one spatial dimenion, the interacting boson with hamiltonian (!69l) is known to 
be exactly solvable using the TBA. In particular, the free energy is 

I (dk) log {l + e-^<'^), (88) 
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where e(/c) obeys the Yang- Yang equation: 

e{k) = k^-^ J {dk')K{k, k') log(l + e-"'^''^), (89) 
The TBA kernel K is the derivative of the log of the S-matrix 

^(M-')^l^i;^f;i (.0) 

K{k,k') = -^d,\ogS{k,k') = j^,^,"/^,^^^^, ^ (91) 



Comparing with (l33l) . one sees that the signs in ( 18811 are fermionic. The TBA is 
really a description of the model in terms of its fermionic dual. 

Expanding the right hand side of flHHl) to the leading order in the kernel K using 
(|89l), we have 



^ ^ J-o + / {dk)fo{k) 



J {dk')K{k, k') log (1 + e-'^'^^') 



(92) 



Integrating the right hand side by parts with respect to k' one obtains 

I {dk){dk')fo{k)fo{k') (^-2z\k - k'\ log (^ 1^ ~ ~ ^^^/^ ^^ (93) 

The above expression coincides with the = 1 term of (jUj) , the first order correction 
of our formulation, with the kernel G2{k, k') given by fl57|) . 
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IX. APPENDIX A 

We shall now show that the saddle point value of F as defined in (l53l) equals the 
free energy JF given in fl37|) . Throughout this section the integral equation fl62|) is 
always assumed to be satisfied. 
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First, let us define 

A/(k)^7(k)-/o(k), (94) 

i.e. A/ is the sum of all digrams with two external legs and any non-trivial insertion. 
In particular, the sum is not restricted to 2PI terms only; all one-particle reducible 
terms are also included. By re-arranging fl56l) . we get one useful identity: 

A/ s{f - /o) 



/o 1 + sfo 

We will now rewrite in a more convenient form using (p5|) : 

1 r.-, / - f f-h 



/3 
1 

1 



(rfk) slog(l + s/) 



sfo 



(95) 



(96) 



(rfk) 



slog(l + sfo) + slog 1 



A/\ A/ 



fo 



— s- 



fo 



^ n=2 



:A/ 

/o 



As for the using interaction part of F, one can show using the integral equation fl62l) : 

1 

'2P 



Fi = 



{dk){dk')f{k)f{k')w2{k,k') 

A/(k) 



(97) 



/o(k) 



Together we conclude that the saddle point value of F is 

A/ 



F 



-^0 + -^ lidk) 



lA/ 
2 /o 



n=2 



fo 



(98) 



We shall now give this equation a graphical meaning. Recall that A/ represents 
the sum of all diagrams with two external legs and some non-trivial structure between 
the two legs. A///o is then the sum of such diagrams with one of the two legs 
amputated. The quantity s J (dk) {Af (k) / foik))"" then represents a ring with n non- 
trivial insertions, with the factor of s again coming from the loop. (See Figure [8]) 
Again we stress that such insertion needs not be IPI. 
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o 



(a) 



(b) 



FIG. 8: (a) Graphical representation of A///o; one of the two external legs is amputated, 
(b) The graph for s J {dk){Af{k)/fo{k))". In this example n = 5. 



The correction to the free energy F — J-'o can then be graphically represented as 
Figure m The factor l/n of the nth order term in fl^S]) is cancelled because there are 
n ways to order these n insertions to form an identical ring, not counting possible 
extra symmetry. 

Recall that, within the foam-diagram approximation, the free energy fl37|) is really 
JFq plus the sum of all connected foam diagrams. Every term in Figure [9] is a foam 
diagram by construction. To prove the equality of ( 137|) and ( 198|) . we only need to 
verify (l98l) . or equivalently Figure [9l produce the correct coefficient for every foam 
diagram. 




+ . . . 



FIG. 9: Graphical representation ol F — J-^o- 
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Consider a foam diagram with loops labeled i = 1 . . . N, and the ith loop has 
rii vertices attached to it. It can be shown that 

N 



2N -2. 



(99) 



i=l 



Now we consider how the particular foam diagram may be constructed from Figure [H] 
if the loops shown is the ith loop of the diagram. The i-th loop has Ui vertices; it can 
be built from terms with at most rij insertions. The m-th order term contributes only 
if m < Ui- The number of ways to partition the IPI insertions into m connected 
groups is simply the binomial coefficient C^. 

We can now sum over every loop of the diagram to get the overall coeficient: 



N 



m=2 



N 

E 

i=l 

N 

E 

i=l 
1 



1 



j2i-ircz--ci+c: 



m=0 



(100) 



We have used the identity (!99|) to evaluate the sum. This shows that the coefficient 
of a diagram without symmetry equals exactly unity. If a diagram has an A^-fold 
rotational symmetry, then the above counting argument overcounts the contribution 
by times, and a symmetry factor 1/A^ is needed. 

We have thus proven that Fi in eq. fl97j) generates every foam diagram with the 
correct coefficient. Its saddle point value is therefore the correct free energy of the 
system within the foam-diagram approximation. 
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